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The nature of the interaction between ultracold atoms with a large orbital and spin angular 
momentum has attracted considerable attention. It was suggested that such interactions can lead 
to the realization of exotic states of highly correlated matter. Here, we report on a theoretical study 
of the competing anisotropic dispersion, magnetic dipole-dipole, and electric quadrupole-quadrupole 
forces between two dysprosium atoms. Each dysprosium atom has an orbital angular momentum 
L = 6 and magnetic moment /i = 10/is. We show that the dispersion coefficients of the ground state 
adiabatic potentials lie between 1865 a.u. and 1890 a.u., creating a non-negligible anisotropy with 
a spread of 25 a.u. and that the electric quadrupole-quadrupole interaction is weak compared to 
the other interactions. We also find that for interatomic separations _R < 50ao both the anisotropic 
dispersion and magnetic dipole-dipole potential are larger than the atomic Zeeman splittings for 
external magnetic fields of order 10 G to 100 G. At these separations spin exchange can occur. We 
finish by describing two scattering models for inelastic spin exchange. A universal scattering theory 
is used to model loss due to the anisotropy in the dispersion and a distorted-wave-Born theory 
is used to model losses from the magnetic dipole-dipole interaction for the ^''^Dy isotope. These 
models find loss rates that are the same order of magnitude as the experimental value. 



I. INTRODUCTION 

In recent years significant effort has been devoted to 
the characterization of the interactions with submerged- 
shell 3d-transition-metal and 4f-rare-earth atoms [iHS]. 
These atoms have an electronic configuration with an 
unfilled inner shell shielded by a closed outer shell. They 
also tend to have a large magnetic moment due to a large 
number of unpaired electrons, which presents opportuni- 
ties to explore the effect of anisotropic magnetic dipole- 
dipole interactions between them. Long-range dipolar in- 
teractions create conditions for realizing novel quantum 
states of highly correlated ultracold atomic matter [71 [9] . 
This physics complements that proposed with ultracold 
polar molecules, another system in which exotic quantum 
phases are predicted [T0Hl3] . Here, dipole-dipole forces 
originate from a non-zero electric dipole moment. Unlike 
for magnetic atoms, however, the electric dipole moment 
must be induced by an external electric field. 

Submerged shell atoms are expected to have sig- 
nificantly suppressed inelastic, energy-releasing spin- 
exchange collisions because of shielding caused by 
the closed outer-shell electrons. This effect was 
first predicted and demonstrated for collisions between 
submerged-shell atoms with helium [JHH HI] • The sup- 
pression of inelastic loss with the He atom indicates that 
there is no coUisional anisotropy. The spherically sym- 
metric He atom can not take up angular momentum from 
the submerged-shell atom. 

Recent measurements of the spin-exchange rates be- 
tween two submerged-shell atoms, however, have seen 
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no suppression and, in fact, the rate coefficients are of 
the same order of magnitude as for non-submerged shell 
atoms O m [15] . A possible explanation for this phenom- 
ena, given in Ref . [5] , is that most submerged-shell atoms 
have a non-zero orbital electron angular momentum L. 
This leads a non-zero electrostatic quadrupole moment 
and anisotropic quadrupole-quadrupole interaction that, 
in principle, can cause substantial losses. 

In this paper we propose and discuss another mecha- 
nism that leads to losses. We will show that the large loss 
rate of order 10~^" cm^/s, observed in [HI |51 [TSj, might 
have been due to anisotropy in the dispersion forces at 
short inter-atomic separations. This anisotropy is also in- 
duced by the nonzero L. We study this new mechanism 
of spin-exchange for the submerged-shell atom with the 
largest magnetic moment, dysprosium. It has an unfilled 
4f^° shell lying beneath a filled 6s^ shell leading to a large 
orbital, L = 6, and total, j = 8, angular momentum. Its 
ground ^Ig state has a magnetic moment oi ^ — lO/is, 
where /is is the Bohr magneton. Only recently, the first 
laser cooling and trapping experiments of a large number 
of dysprosium atoms have been reported jl5j. The first 
measurements of inelastic coUisional rates in this study 
suggest that anisotropy in the inter-atomic forces plays 
a significant role. 

The paper is organized as follows. We first analyze 
the isotropic and anisotropic dispersion interaction be- 
tween two Dy atoms in Section |ll] and compare it with 
the magnetic dipole-dipole and electrostatic quadrupole- 
quadrupole interactions. The dispersion coefficients are 
calculated from atomic transition frequencies and dipole 
moments. The quadrupole moment of Dy is determined 
from a multi-configuration electronic structure calcula- 
tion. In Section |III| we study the relative strength of the 
interactions in the presence of an external magnetic field 
and rotation. In Section ITVl we use these interactions to 
find the first estimates of the spin-exchange loss rates and 
compare with experimental results. 
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II. RELATIVE STRENGTH OF INTERACTION 
FORCES BETWEEN GROUND STATE DY 
ATOMS 

The theoretical calculation of the ground state Dy2 
potentials and their dispersion parameters is challeng- 
ing due to the complexity of the spin structure of the 
ground-state ^Ig Dy atom. For example, there are 81 
gerade and 72 ungerade potentials that dissociate to the 
^Ig -I- ^Ig limit. In spite of this complexity we have be- 
gun to calculate the van der Waals Cg coefficients for two 
interacting Dy atoms. 

For two colliding atoms we can define the angular mo- 
mentum J = Ji -l-j2, its projection M along the direction 
of the external magnetic field B, and its projection 
along the internuclear axis. For this relativistic molecule 
the adiabatic Born-Oppenheimer (BO) potentials are la- 
beled by n^, where a = g/u for gerade and ungerade 
states, respectively. Gerade (ungerade) symmetry corre- 
sponds to superpositions of even (odd) values of J. The 
superscript ± is only relevant for = states. For each 

there are 17 — \Q\ adiabatic Born-Oppenheimer (BO) 
potentials combined. 



where eq is the electric constant, di and d2 are the electric 
dipole operators for the two atoms, and diz and d2z are 
their components along the internuclear axis. 

Using the Wigner-Eckart theorem we write the matrix 
Cg as 
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A. Electrostatic dispersion interaction 

We describe the dispersion interaction potential for two 
ground-state atoms in state |jimi,j2m2) using degen- 
erate second-order perturbation theory similar to that 
given by Ref. |16|. The magnetic quantum numbers mi 
and m2 are projections along the internuclear axis of 
the total atomic angular momenta ji and j2 for the two 
atoms, respectively. Here ji = 22 = 8. (In this section 
we break with convention and use roman symbols for 
atomic projection quantum numbers on the internuclear 
axis.) Matrix elements of the dispersion potentials are 



(jl™l,i2"^2|^7disp|Jln^^, J2TO2) 
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(jin^l , J2'7l2 I Vdd\naj atria , jfcmfc) 

X {uaja rria , Ubjbmb I Vdd \ji m'l , j2 1^2) , 

where the Ce{mim2,mim2) form a matrix of dispersion 
coefficients, R is the separation between the atoms, the 
sums are over all electronic states \naja'm'a,n'bjbmb) of 
atoms a and b excluding states with energies En^j^ and 
Enbjb equal to the ground state energies Ei and E2 ■ The 
operator Vdd is the dipole-dipole interaction Hamiltonian 
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Note that the A-'^-'^-'"-"' , , conserve the molecular vyvo- 

mim2,m^m2 ^ 

jection = mi -|- m2 = m'l -|- m2 and are independent on 
atomic transition frequencies and dipole moments. For 
this homonuclear molecule gerade/ungerade symmetry 
states are most conveniently constructed by transform- 
ing to states of total J. That is to states \{j\j2)J^ and 
noting that even(odd) J states have gerade(ungerade) 
symmetry. 

There are six independent i^j^j^ for two Dy ^Ig atoms 
as the selection rules of the electric dipole operator re- 
quires that |ii-l| < ja < Ji + 1 and |j2-l| < jb < j2 + l- 
For homonuclear dimers the ^^j^j^ is symmetric under 

interchange of ia and jb- We have determined K-'^^^^ us- 
ing 62 experimental transition frequencies and oscillator 
strengths from the ground to various excited states of the 
Dy atom [T7]. Table |l] lists the values of Kf^^j^- 

The adiabatic dispersion potentials and, thus, the long- 
range of the Born-Oppenheimer potentials are found by 
the diagonalizing the Cg matrices for each Figure [l] 

shows the adiabatic gerade and ungerade Cg coefficients 
as a function of the projection quantum number of the 
total angular momentum J on the interatomic axis. The 
number of of adiabatic Cg values is smaller for larger 
57. In fact, for il = 16 there is only one potential. It 
has gerade symmetry. In total there are 81/72 disper- 
sion coefficients corresponding to the ground state ger- 
ade/ungerade potentials. The coefficients in Fig. [l] show 
a smooth nearly parabolic behavior with the projection 
number fi. This dependence is a consequence of the 
anisotropic coupling of the open /-shell electrons of the 
two atoms. As a result the interaction energy depends 
on the relative orientation of the atoms. 
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FIG. 1: Gerade (filled circles) and ungerade (open circles) adi- 
abatic Ce coefficients for the interaction between two ground 
^Ij=8 state Dy atoms as a function of the projection Q of 
the total angular momentum J on the interatomic axis. The 
difference between the dispersion coefficients for the ger- 
ade/ungerade symmetry is small and invisible on the graph. 
A larger Ce coefficient implies a deeper Born-Oppenheimer 
potential. 



FIG. 2: Gerade (filled circles) and ungerade (open circles) adi- 
abatic C3 coefficients for the interaction between two ground 
^Ij=8 state Dy atoms as a function of the projection Q of the 
total angular momentum J on the interatomic axis. A larger 
C3 coefficient implies a deeper Born-Oppenheimer potential. 



B. Magnetic dipole-dipole and electrostatic 
quadrupole-quadrupole interaction 

The matrix elements of the magnetic dipole-dipole in- 
teraction between two magnetic dipoles /Z = gjHsj is 



{jimi , J2m2 1 C/mdd \jim[ , ^2^2) 



C3{mim2,m[m2) 

i?3 



= {jimi, j2m2\Vf,f,\jim[, j2m2) , (4) 



where Vnu is magnetic dipole-dipole operator 



An 



i?3 
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and gj = 1.24159 is the g-factor for the ground ^Ig state 
of the Dy atom [T8|, and /io is the magnetic constant. A 
more accurate value for the magnetic moment of Dy is fi 
= gjiJ'B X j = 9.93 ^s- 

Figure [2] shows the adiabatic gerade and ungerade C3 
coefficients as a function of SI. These coefhcients are ob- 
tained by diagonalizing the matrix Eq. Q at each R. The 



TABLE I: The J^'j^j^ matrix elements in atomic units for the 
dipole transitions from j\ — j2 ~ 8 to ja, jb ~ 7, 8, or 9 for 



two interacting Dy atoms. 
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FIG. 3: Adiabatic gerade interaction potentials of the com- 
bined electrostatic dispersion and magnetic dipole-dipole 
forces between two Dy atoms in the ground ^Ig state and 
projection Q=0. Here ks is the Boltzmann constant. The 
effect of rotation is not included. 



values are both positive and negative. A comparison with 
the adiabatic Cg coefficients in Fig. [T] shows a different H. 
dependence. 

A more accurate description of the long-range interac- 
tion is obtained by first adding the dispersion t/disp and 
magnetic dipole-dipole [/mdd interactions together and 
diagonalize at each internuclear separation R. Unlike, 
for the previous cases the eigenfunctions now depend on 
R. As an example, the resulting adiabatic gerade poten- 
tials for projection Q—0 as a function of R are shown in 
Fig. [3] At small R the dispersion interaction dominates. 
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whereas for R > 150 aq the magnetic dipole-dipole in- 
teraction plays a major role. For intermediate R these 
forces compete leading to both attractive and repulsive 
potentials depending on sign of the C3 coefficient. 

Our unrestricted coupled cluster calculation with 
single, double, and perturbative triple excitations 
UCCSD(T) [T5] shows that the quadrupole moment of 
the Dy atom in the ^Ig state is very small and equal to 
(5^-0.00524 a.u.. As a result the quadrupole-quadrupole 
interaction energy is seven orders of magnitude weaker 
than the other atom-atom interactions. 



III. INTERACTIONS IN A MAGNETIC TRAP 

We now analyze the relative strength of all interac- 
tions between two Dy atoms in a magnetic field. The 
magnetic field is added as either the atoms are held in a 
magnetic trap ^5j with a spatially varying field strength 
or are held in an optical trap with an homogeneous B 
field to control the interaction between the atoms. In 
addition, the molecule can rotate, which is described by 
the Hamiltonian U^i^ /{2mrR'^), where £ is the relative 
orbital angular momentum between the two atoms and 
rrir is the reduced mass. 

It is convenient to choose a coordinate system with pro- 
jection quantum numbers defined along the external mag- 
netic field direction. Again following convention, projec- 
tion quantum numbers are labeled by roman symbols. In 
this coordinate system the rotational and Zeeman inter- 
actions as well as the isotropic or "average" dispersion 
potential shift molecular levels, whereas the magnetic 
dipole-dipole interaction and anisotropic component of 
the dispersion potential lead to coupling between differ- 
ent rotational and Zeeman components. As a result, the 
angular momentum projection M of J can change up to 
2 units due to the magnetic dipole interaction and up to 
4 units due to the anisotropic dispersion potential PU] . 

Figure [4] shows various anisotropic properties that can 
lead to reorientation of the Dy angular momenta as a 
function of R. Firstly, the Zeeman splitting gj^gB be- 
tween neighboring magnetic sublevels for magnetic field 
strength of 10 and 100 Gauss are shown. The anisotropic 
potential ACq/R^ is drawn assuming a typical value 
of ACe = 25 a.u., based on the spread of Cq value 
shown in Fig. [T] We also present the splitting between 
the rotational levels ^=0 and 2 of the ground state as 
/{2mrR?)- Finally, the value of splitting due to mag- 
netic dipole-dipole interaction is given. 

For different interatomic separations different forces 
dominate. In fact, when the curves for the magnetic 
dipole or anisotropic dispersion interaction cross the Zee- 
man or rotational energies spin flips can occur. At large 
R the Zeeman splitting is dominates. Both magnetic 
dipole-dipole and anisotropic electrostatic curves cross 
the Zeeman B =100 G curve at i? < 35ao, where chem- 
ical bonding should play an important role as well. For 
the weaker magnetic field of i? = 10 G the spin coupling 
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FIG. 4: Level splitting due to the dominant interaction forces 
in atomic units as a function of interatomic separation. In 
atomic units the Zeeman splitting is gjB/2, the splitting be- 
tween (.—Q and 2 rotational levels is given by &/{2mrR^), the 
splitting due to the magnetic dipole-dipole (MDD) interaction 
is 2oi^ j{gj /2)^ / E? , where a is the fine structure constant, and 
the anisotropic dispersion (AD) interaction is ACe/R^, where 
AC6=25 a.u.. Here rrir is the reduced mass in units of the 
electron mass. 




FIG. 5: Dispersion interaction potentials for Cs ~ 1890 a.u. 
for the lowest five partial waves with centrifugal barriers for 
p-, d-, /-, and g- waves 

occurs for R near 50ao. The interactions will lead to 
mixing of rotational levels for R < 50ao as well. 



IV. FIRST ESTIMATE OF INELASTIC RATE 
COEFFICIENTS 

In this paper we complete our analysis of the interac- 
tion between ultracold ^^^Dy atoms by a first estimate 
of inelastic loss rates due to the anisotropy in the dis- 
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persion and magnetic dipole-dipole interaction. We will 
perform a separate estimate of losses from these interac- 
tions. For both cases we consider atoms in a magnetic 
field with a strength of the order of 1 G as described in 
the recent experiment [T5^, with the goal to model its 
losses. The experiment started from a gas of Dy atoms 
in a quadrupole magnetic trap with atoms distributed 
over atomic magnetic sublevels with positive magnetic 
moment. le. states \jm) with m > where the pro- 
jection m is defined along the magnetic field direction. 
Inelastic spin-exchange collisions to states with m < 
lead to atom loss. 

We first describe a model for the loss due to the 
anisotropic dispersion potentials based on an universal 
single-channel scattering model developed and used in 
Refs. pT| - B5] . This universal loss model assumes scat- 
tering from a single potential of the form —Cq/R^ + 
h^l{i + l)/(2TO,.i?2) for R > Rc and that all flux that 
reaches the critical separation Rc undergoes irreversible 
spin-exchange independent of collision energy and partial 
wave i. 

For Dy we can use this universal model under several 
assumptions. We first note that the anisotropy ACg of 
the dispersion potential is small compared to the average 
or isotropic dispersion potential. Secondly, an external 
magnetic field is applied that splits the different m levels 
and, as shown in Fig. [4j the spin flip occurs between 
Rc = 35ao and 50ao depending on the magnetic field 
strength. We can therefore apply the universal model 
assuming a mean isotropic Cg value and that, due to 
the anisotropic dispersion potential, no flux returns from 
R < Rc- 

For temperatures between 100 fiK to 1 mK only a few 
partial waves £ contribute to the collisions. Figure [5] illus- 
trates this by showing the centrifugal barriers for s,p^d^ f 
and g partial waves as a function of R. The tempera- 
ture range of interest lies well below the g-wave barrier. 
Within the universal scattering model the contribution 
to the inelastic rate coefficient for partial wave i and 
projection is 

Kim^E) ^V~{1~ \SirnAE)\'') , (6) 

where E — kf/{2mr) is the collision energy, ki is the 
initial relative wavenumber, Vi is the initial relative ve- 
locity, and the Simg (E) are diagonal scattering S'-matrix 
elements. The solution 'i'emi (R) of the radial Schrodinger 
equation for the single-channel potential with boundary 
condition 

at short range R < Rc with R^ ~ \j imrC^/h? and 

p-ikiR ^ikiR 
^hne {R) = 7= Slmi {E) 

at large R determines the Sirai matrix elements. The 
partial and total loss rate coefficient are f3e{E) = 




E/k„ (nK) 



FIG. 6: The inelastic loss rate coefficient for a non-spin- 
polarized sample of ground state ^^*Dy atoms as a function 
of collision energy based on an universal scattering model 
for losses due to the anisotropy of the dispersion poten- 
tial and a Born approximation for losses from the magnetic 
dipole-dipole interaction. For the universal scattering model 
rate coefficients for the lowest four partial waves as well as 
summed rate are shown. The loss rate coefficient for the mag- 
netic dipole-dipole interaction is indicated by the abbreviation 
MDD. The unitary limited loss rate coefficients for the lowest 
four partial waves are plotted as dashed lines. 



^Em.Kim.iE) and (3{E) = 2 (i?), respec- 
lively. The factor of 2 is due to the fact that after each 
collision two atoms are lost. Partial and total inelastic 
loss rate coefficients are shown in Fig. [6] for collision ener- 
gies upto 1.5 mK. The figure shows that the loss rate for 
the different partial waves becomes large for collision en- 
ergies approaching the corresponding centrifugal barrier. 
Moreover, except for extremely small collision energies 
the total loss rate coefficient slowly increases with en- 
ergy. For comparison we have also indicated the unitar- 
ity limit for each partial wave. The unitary limit occurs 
when \Sernf{E)\'^ = for all collision energies. 

We now turn to a model for losses due to spin exchange 
induced by the magnetic dipole-dipole interaction. For 
simplicity we assume that the atoms are in the stretched 
state with m — +j. We estimate this inelastic loss rates 
using first-order perturbation theory similar to that ap- 
plied for the calculation of the dipolar relaxation rates in 
a gas of chromium atoms [T]. We immediately note that 
the magnetic moment of dysprosium or chromium atoms 
is large and, therefore, a perturbative theory may not 
provide an accurate loss rates. However, it is expected 
to give a reasonable estimate. Following Ref. [1 , the loss 
rate coefficient for a single spin flip with M — ^ M — 1 
(i.e. from M=16 to 15), averaged over all possible rela- 
tive orientations of the initial relative momentum fc^, is 

i.-Y^j [ ) [i + MV^O]— , (7) 
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where h'^k'j/(2fir) = gjl^sB and 

hix) - -l/2-(3/4)(l-x2)2log[(x-l)/(x+l)]/(x(l+x2)) 
for X > 1. 

Similarly, for a double spin flip with M — > Af — 2 the 
rate coefiicient is 

72^j^J [ ^^^2 ) [l + h{k,/h)] — , (8) 

where now fi^k^/ {2fj,r) = 2gjfXBB. The total dipole- 
dipole loss rate is given by 7 = 2(71 + 72) and shown 
as a function of collision energy in Fig. |6] 

The loss rates in Fig. |6] are smaller than the rate mea- 
sured in Ref. [IJl as 2.1(2) x 10~^° cnr" /s for temper- 
atures around 500 /iK. This suggests the presence of a 
resonance in the scattering process. In fact, flux of atoms 
can return from small R, interfere with the incoming flux 
to lead to an increasing loss. We obtained a similar ef- 
fect in our analysis of the reactive collisions between two 
KRb molecules [21]. 

V. CONCLUSION 

In conclusion, we have studied the origin of the 
anisotropy in the long-range interaction between ground 
state dysprosium atoms. This is a first step towards a 



complete multi-channel description of inelastic and elas- 
tic collision between such atoms. We find van der Waals 
coefhcients by using known atomic dipole moments and 
energy levels. Our coefficients form a lower bound. We 
show that the splitting between or anisotropy of the 
Born-Oppenheimer potentials is almost two order of mag- 
nitude smaller than their average or isotropic potential. 
In addition, we have presented two approximate single- 
channel calculations to estimate inelastic losses when 
Dy atoms are not in the energetically-lowest Zeeman 
sublevel. The first model describes losses due to the 
anisotropy of the dispersion potentials and is based on 
an universal scattering theory. The second perturbative 
model describes losses due to the magnetic dipole-dipole 
interaction. The only way to obtain a clear and quantita- 
tive understanding of collisions between Dy atoms is by 
a coupled-channel calculation. We will do so in the near 
future. It will enable us to predict location of magnetic 
Feshbach resonances in the energetically-lowest Zeeman 
level. 
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